Equilibrium Correlation Functions of the Spin-Boson Model with Sub-Ohmic Bath 
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The spin-boson model is studied by means of flow equations for Hamiltonians. Our truncation 
scheme includes all coupling terms which are linear in the bosonic operators. Starting with the 
canonical generator r]c = [H(i,H] with H^ resembling the non-interacting bosonic bath, the flow 
equations exhibit a universal attractor for the Hamiltonian flow. This allows to calculate equilibrium 
correlation functions for super-Ohmic, Ohmic and sub-Ohmic baths within a uniform framework 
including flnite bias. 
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The spin-boson model, consisting of a spin-l/2-system 
coupled to a bath of non-interacting harmonic osciUa- 
tors, is the simplest non-trivial model that captures the 
essential influence of dissipation on a quantum mechan- 
ical system 0, |^. The physical behavior of the system 
is determined by the low-energy modes of the bosonic 
bath. Modeling the interaction such that one recovers 
the classical phenomenological equation of motion with 
a frictional force proportional to the velocity, one speaks 
of Ohmic dissipation. At zero temperature the system 
then undergoes a continuous transition from delocaliza- 
tion to localization at a critical interaction strength. If 
the low-energy modes are suppressed or enhanced rel- 
ative to the Ohmic bath one speaks of super-Ohmic or 
sub-Ohmic dissipation. For super-Ohmic dissipation, the 
effect of the bath on the system is generally that of renor- 
malizing the system parameters. For sub-Ohmic baths, 
the system has the tendency to localize. 

Physical systems mostly exhibit super-Ohmic or 
Ohmic dissipation. Obvious realizations for sub-Ohmic 
baths are not known and the only known results are valid 
for high-temperatures and small time scales They 
are e.g obtained using Monte Carlo simulations |^ . The 
lack of results for long time scales at zero temperature is 
mainly due to technical difficulties since the sub-Ohmic 
coupling represents a relevant perturbation in the lan- 
guage of the renormalization group . 

With the upcoming interest in quantum computing 
based on solid state qubits, investigations on different 
types of dissipative systems have emerged ^, 0] and 
there are e.g suggestions to model 1// noise found in re- 
cent experiments by a limiting sub-Ohmic bath 
The widespread believe that there is total suppression 
of quantum tunneling at zero temperature for any sub- 
Ohmic bath independent of the coupling strength would 
pose a serious problem to quantum computing within a 
solid-state matrix due to 1// noise. But results based 
on the flow equation approach pO| and on exact math- 
ematical techniques indicate that there is a first or- 
der phase-transition from delocalization to localization at 
some finite coupling constant. To our knowledge, spec- 



tral properties which are valid on all energy scales have 
not been presented yet. 

In this letter, we will investigate the spin-boson model 
by means of flow equations since this approach is 
not confined to any particular bath type. We are thus 
able to treat all environments on the same footing. One 
crucial difference to prior works on the method applied 
to the spin-boson model is, that we allow all coupling 
terms which are linear in the bosonic bath operators to 
be generated. Furthermore, we set up flow equations 
that exhibit universal asymptotic behavior which assures 
a systematic decoupling procedure where the low-energy 
modes are decoupled last. In view of the manipulation 
of qubits, we also treat the spin-boson model with finite 
bias. This gives rise to distinguished bosonic modes and 
it will also be necessary to vary the unitary representation 
of the initial Hamiltonian. 

We start our discussion with the reflection-symmetric 
model where the initial Hamiltonian is given by 
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Notice that summation over the bath modes a is implied 
throughout this work. The objective is to evaluate the 
equilibrium correlation function C{t) — {{ijx{t),ax\) /2. 

The scale of the isolated system shall be given by 
Aq = A(£ = 0). The coupHng constants are as usual 
parametrized in terms of the spectral coupling function 
J{uj) = {\lf5{uj - (jJa) = 2aifi-«w"/(w/wc) with the 
coupling strength a, the bath cutoff ujc and the cutoff 
function f{x). The different bath types are character- 
ized by s, where s = 1 stands for Ohmic baths and s > 1 
and s < 1 for super- and sub-Ohmic baths, respectively. 
Non-Ohmic baths also give rise to an additional energy 
scale K. 

The Hamiltonian is now subjected to a series of in- 
finitesimal transformations which can be written in differ- 
ential form as diH = [77, H], £ denoting the flow parame- 
ter and resembling the square of the inverse energy scale 
which is being decoupled The anti-Hermitian gen- 
erator 77 is first chosen to be canonical, i.e. 77c — [Hq, H] 



|r2[ , with Hq = Uabl^ba. This leads to the trivial fixed 
point Hamiltonian H(£ = oo) = Hq. But r]c generates 
a new coupling term which is linear in the bosonic oper- 
ators and given by Hnew — io'y^{ba ~ ^a)- This term 
is included in the fiow of the Hamiltonian and also con- 
tributes to the canonical generator. 

Moreover, the commutator [77^ H] contains coupling 
terms which are bilinear in the bosonic operators. These 
terms are not included in the Hamiltonian flow explicitly, 
but we truncate the Hamiltonian after coupling terms 
linear in the bosonic operators. Yet, we consider bilin- 
ear coupling terms adequately by introducing a new term 
to the generator. Doing so, we are guided by the exact 
solution of the dissipative harmonic oscillator via flow 
equations p3[ |. The additional term of the generator is 
given by ijnew = ila,a'{ba + bl,){ba' -b^,), where the pa- 
rameters ?7q,q' are to be chosen such that the bilinear 
terms are not being generated. This cancellation cannot 
be done exactly, but we will neglect operators which are 
normal ordered with respect to the ground-state of the 
one-particle Hamiltonian HP — and the bath Hamil- 
tonian Hq, respectively. The parameters rja^a' for a ^ a' 
are then given by (jia,a = 0) 
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The renormalization of the bath modes cua vanishes in 
the thermodynamic limit and will be neglected. For 
the tunnel-matrix element we obtain the flow equation 
dgA = —2uJa^%^a- The flow of the couphng constants 
is described by the following non-linear differential equa- 
tions: 



(2) 



In order to calculate correlation functions the observ- 
able has to be subjected to the same sequence of infinites- 
imal transformations as the Hamiltonian. The truncated 
flow of the observable shall be given by ax{t) — h^{£)ax + 
(^zXa{^){ba + bl). The flow equations dga^ = [ri,(Tx\ are 
closed according to the same normal ordering scheme as 
above and thus read 



Characteristic for dissipative systems is that the one- 
particle parameter ^ for £ 00, i.e. the spec- 
tral weight of the system dissipates into the bath [ p^ . 
The correlation function C{t) = dcoe^'^^C^uj) is then 
determined by the spectral function C{uj) = {Xai^ — 
oo))^5(w — uja) which is obtained by numerical integra- 
tion of the flow equations. In Fig. |l|, C{uj) is shown for 
various bath type s. The parameters of the spectral cou- 
pling function J{uj) are chosen as a = 0.1, Wc/Aq ~ 10, 
K/Aq = 1. The spectral functions yield the correct low- 
energy behavior C(w) oc for w/Aq <C 1 which clearly 
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FIG. 1: The spectral function C{uj) for J{i^) = 
2aK^-"uj'e{ujc - Lo) with a = 0.1, lJc/Aq = 10, K/Aq = 1, 
and e = for various bath type s. 



shows the non-perturbative nature of our approach. We 
will come back to this point later. 

For the numerical implementation of the flow equa- 
tions, we employed a linear energy spacing, which allows 
for an easy evaluation of the principal value integral aris- 
ing from rja,a' of Eq. (|l|), and the sharp cutoff function 
f{x) — Q{1 — x). The results for low energies do not 
depend on the specific choice of f{x) nor on the energy 
discretization. For larger tOc a logarithmic energy spacing 
is needed which yields equivalent results. 

In Fig. ^, the normalized spectral function S{uj) cx 
C{ijj)/lj'^ is shown for various coupling strength a where 
cj = A* denotes the maximum of C(w), i.e. A* « 
Ao(Ao/ti-'c)"^'-^^"'' for Ohmic baths. The parameters of 
J([x)) are chosen as ujc/ Aq = 10, K/Aq = 1, s = 1 (left 
hand side) and s — 0.8 (right hand side), respectively. 
The transition from coherent to incoherent tunneling oc- 
curs at the coupling strength a* where the resonance at 
Lo = A* disappears. We obtain q;*>0.4 for an Ohmic 
bath and a* w 0.3 for a sub-Ohmic bath with s = 0.8. 

Our truncation scheme is valid for small and interme- 
diate coupling strengths a<0.4, and for Ohmic baths our 
approach yields results similar to those obtained by Costi 
and Kieffer by means of the numerical renormalization 
group (NRG) |l4j. But since NRG is based on the map- 
ping of the spin-boson model onto the anisotropic Kondo 
model it cannot be used for sub-Ohmic baths. 

We will discuss our approach in more detail, now. 
We have chosen the diagonal Hamiltonian to be Hq = 
^cMyba- This choice indicates that A(f) ^ as £ ^ 00, 
i.e. the fixed point Hamiltonian is scale-independent and 
thus universal. With no asymptotic scale present, the 
flow equations will systematically first decouple high- 
energy modes and finally low-energy modes. This de- 
coupling procedure is passed on to the observable flow. 
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FIG. 2: The spectral function S{u}) for J{uj) = 
2aK^-"Lo''e{ujc - tj) with cj^/Ao = 10, K/Ao = 1, e = 0, 
Ohmic bath s — 1 (left hand side) and sub-Ohmic bath 
s = 0.8 (right hand side) for various coupling strength a. 



i.e. the spectral function C(lu) is first determined at high 
energies and the trivial low- frequency behavior of the sys- 
tem is calculated at last. Moreover, the spectral function 
around the resonance uj w A* is determined by the stable 
flow of intermediate i « (A*)~^, away from the asymp- 
totic regime of decoupling quasi-degenerate states. An- 
other advantage of letting the tunnel-matrix element A 
rcnormalize to zero is that the fixed point Hamiltonian 
Hq possesses the correct spectrum, consisting of a non- 
degenerate ground-state and a gapless continuum. 

The flow equations can numerically be integrated only 
up to i* « {5uj)~'^ with Sui denoting the linear en- 
ergy spacing. Determining the asymptotic behavior 
will thus specify the spectral function C(w) in the low- 
energy regime which is not accessible by numerical in- 
tegration. A first approximation for the asymptotic 
flow is given by « A^, which is consistent with 
. Introducing the ^-dependent spectral coupling func- 
tion J{uj,£) = (A^(^))^(5(cj — oja), the asymptotic flow 
equations are solved by the ansatz A{£) — > a£^^^'^ and 
jIuj,£) e-'^^^J{ujV£)- This leads to the following non- 
local differential equation for J{y) with y = uj\/£: 



dyJiy)^-4.iy-a)J{y) + 8J{y) 



ay 7 + t , 



y-y 
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where t = 1 — 4 dy' J{y') and the boundary condition 
« = dyyJ{y). 

The asymptotic behavior of the observable flow is de- 
scribed by a similar ansatz, i.e. h^{£) M^^^^^'^ 1^ 
and S{u:,£) = XyMxUm^-^c) ^ £-^/^-'' /^Siu;V£). 
The non-local differential equation for S{y) is similar to 
the one for J{y) |l^. Defining the ^-dependent spectral 
function C(u;,i!) = {Xai£))'^S{uj—uja), we obtain for be- 
ing in the asymptotic regime C{uj) — C{uj,£*) = lo^ I{y*) 



with y* = Loy/I* and 
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For an initial spectral coupling function J{lo) cx o;*, the 
numerical analysis yields s ~ s' independent of the cou- 
pling strength. The asymptotic scale-invariant flow thus 
induces the correct long-time behavior of the correlation 
function C{uj). 

There are two criteria for the assessment of the quality 
of the results. One is the normalization of the spectral 
function C{uj), i.e. diijC{uj) — 1. The second is the 
Shiba relation lim^^QC{uj)/uj'^ — 2aK^~^ C{uj)/uj 
|]l5| . The numerical results show that the spectral sum 
rule is fulflUed within less than 0.1 % relative error. The 
interpretation of the Shiba relation is not as obvious. 
From (H) we formally obtain lim^^o C{u!)/u!^ — /(O), 
which is not defined since the term dyS{y)/y^~^" di- 
verges for y* ^ 0. But this is not the correct result 
because the proper order how to perform the limites is 
first ^ — !■ 00 and then w — s- 0. The fact that the asymp- 
totic and the thermodynamic limit do, in general, not 
commute is well-known 
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Performing the limits correctly, the left and right hand 
side of the Shiba relation can be determined. This yields 
satisfactory results with relative errors of around 10% for 
Q!<0.2 and s>0.8. For larger a and lower s the relation 
still holds with a maximal relative error of around 30%. 
Approximation schemes like NIBA or numerical methods 
based on Monte Carlo cannot be used to verify the Shiba 
relation since they fail to predict the correct long-time 
behavior. Results from NRG |14 for the Ohmic bath 
have a maximal relative error of around 10%. 

Applying an external bias, the Hamiltonian flow as well 
as the operator flow is not confined to one or two basic 
interaction terms anymore. In this case, we allow all 
coupling terms which are linear in the bosonic operators 
with real coefficients to be generated. The full truncated 
Hamiltonian thus reads 
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with Xa = (6q -I- b\) and pa = —i(ba — bl^). Notice that 
summation over j — e,x, z with cTg = 1 is implied from 
now on. 

In the above representation, the considered equilibrium 
correlation function is defined by C(t) = {{az{t), az})/2. 
The extended truncated observable flow is thus given by 

a, {£) = {£)a, + a^x^, - ayxl ii)Po. , 

where Xai^) = for all £. It is important that the 
constant term h'^ is generated which indicates reflection- 
symmetry breaking, i.e. (ctz) = h'^{£ = 00) = h^. Hence, 
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FIG. 3: The spectral function C{u) for J(<^) = 
2aK^-'' uj" Q{ujc - i^) with a = 0.1, lJc/Aq = 10, A'/Ao = 1, 
and sub-Ohmic bath s = 0.8 for various bias e. At = 
there are (5-peaks with weig ht {Kf (Indent). 



it contributes to the spectral function C{ijj) via a 5- 
function at cj = with weight {h%)'^ . 

In principle, each unitary equivalent initial Hamil- 
tonian would yield the same physical result. Matters 
change as soon as approximations are involved since the 
neglected terms become differently important depending 
on the representation of the initial Hamiltonian. Inves- 
tigating the system with finite bias e we thus introduce 
another parameter 9 which indicates an initial shift of 
the bosonic operators, i.e be ^ be + ^2Zf~' choose 
9 such that the sum rule {crl{£)) = 1 is fulfilled best for 
all £. We also diagonalize the initial Hamiltonian with 
respect to the two-dimensional Hilbert space before we 
apply the flow equations. 

To set up the flow equations for the extended model, we 
proceed in the same way as in the case of the symmetric 
one. Only, the normal ordering procedure now has to be 
defined with respect to distinguished, ^-dependent modes 

ba = ba + 5a/'2' with 

where the one-particle expectation values are evaluated 
with respect to the effective one-particle Hamiltonian 



In Fig. ^, the spectral function C(a;) for initial cou- 
pling J{uj) with a — 0.1, Wc/Aq = 10, K/Aq = 1, and 
sub-Ohmic bath s — 0.8 is shown for various bias e. 
The initial shift of the bosonic operators is taken to be 
9 = e/Ao- This yields a spontaneous breaking of the re- 
flection symmetry which is favorable for the system flow 
as well as for the observable flow. In the indent, the ex- 
pectation value (az) = hi is shown. Similar results are 
obtained for super-Ohmic and Ohmic baths. 

In summary, we have established a flow equation 
scheme which yields reliable results for a large regime 
of the parameter space of the spin-boson model with es- 
sentially new results for sub-Ohmic baths and flnite bias. 
Moreover, we are not conflned to specific cutoff functions 
which might be crucial in view of environments based on 
a structured bath . 
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by the Deutsche Forschungsgemeinschaft. 



HP = -^a, + ^a, with A' = A + A^A^/t^c - KK/^a 
and e' = e — A^A^/wq, — A||jA^/wQ. The distinguished 
modes are derived by quadratic completion and suitable 
factorization, see Ref. [0 for details and extensions. 
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